Quantum spin metal state on a decorated honeycomb lattice 
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We present a modification of exactly solvable spin-| Kitaev model on the decorated honeycomb 
lattice, with a ground state of "spin metal" type. The model is diagonalized in terms of Majorana 
fermions; the latter form a 2D gapless state with a Fermi-circle those size depends on the ratio of 
exchange couplings. Low-temperature heat capacity C(T) and dynamic spin susceptibility x(w,T) 
are calculated in the case of small Fermi-circle. Whereas C(T) ~ T at low temperatures as it is 
expected for a Fermi-liquid, spin excitations are gapful and x(oj, T) demonstrate unusual behaviour 
with a power-law peak near the resonance frequency. The corresponding exponent as well as the 
peak shape are calculated. 
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Quantum spin liquids [lH5| present very interesting ex- 
amples of strongly correlated phases of matter which do 
not follow the classical Landau route: no local order pa- 
rameter is formed while the entropy vanishes at zero tem- 
perature. Although quite a number of different proposals 
for the realization of spin liquid states are available, and 
many interesting results were obtained numerically, see 
e.g. [J], the progress in analitical theory was hindered for 
a long time due to the absense of an appropriate quan- 
tum spin model exactly solvable in more than one spa- 
tial dimension. Seminal results due to A. Kitaev Q may 
pave the way to fill this gap. Kitaev proposed spin-i 
model with anisotropic nearest-neighbour spin interac- 
tions J a afo-" (where a = x, y, z) on the honeycomb lat- 
tice. The model allows for the two-dimensional general- 
ization of the Jordan- Wigner spin-to-fermion transforma- 
tion and can be thus exactly diagonalized. Kitaev model 
realizes different phases as its ground-states, the most in- 
teresting of them corresponds to the region around sym- 
metric point J x = J y = J z = J' . The excitations above 
this ground-state constitute a single branch of massless 
Dirac fermions (Kitaev presents them in terms of Ma- 
jorana fermions those spectrum contains two symmet- 
ric conical points). Thus the Kitaev honeycomb model 
presents exactly solvable case of the critical QSL. Al- 
though long-range spin correlations vanish in this model 
presisely due to its integrability 0], it presents valuable 
starting point for the construction of analitically control- 
lable theories possessing long-range spin correlations. 

Still, here we discuss another extension of the 
model Q: it is interesting to find exact realizations of 
other types of QSLs, the one with gapful excitation spec- 
trum, and the one with a whole surface of gapless exci- 
tations [8 
Kivclson 



Gapful QSL was recently found by Yao and 
They proposed specific generalization of the 
Kitaev model, where each site of the honeycomb lattice 
is replaced by a triangle (we will call it, for brevity, "3- 
12 lattice"), with internal coupling strengths equal to J 
and inter-triangle couplings equal to J'. Topologically 
equivalent structure of such a lattice is shown in Fig. [T] 
Yao-Kivelson model is exactly solvable and contains a 
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FIG. 1: (Color online) "Brick- wall" representation of the dec- 
orated honeycomb lattice. Elementary cell for the Majorana 
Hamiltonian contains 12 sites and is bounded by the blue line. 



critical point at g = J' /J = s/3. At g = y/3 the ex- 
citation spectrum has single low-energy branch of Dirac 
fermions, like the original Kitaev model, whereas at any 
other g the excitation spectrum is gapful. Yao and Kivcl- 
son have shown that the ground state at g < \/3 is a topo- 
logically nontrivial chiral spin liquid with Chcrn number 
C = ±1, whereas g > V3 phase is topologically trivial, 
with C = 0. Exactly solvable QSL model of spin-metal 
type with spins- 1 was proposed by Yao, Zhang and Kivcl- 
son [h]]. In the present Letter we take a different route: 
we show that slight modification of the Hamiltonian of 
spin-i Yao-Kivelson model on the 3-12 lattice leads to 
the spin-metal QSL with a pscudo- Fermi-circle. Similar 
approach was proposed recently in Rcf. 



Ill where QSL 



with pseudo- Fermi-line was found in the spin-^ model 
on a decorated square lattice. Apart from another lat- 
tice studied, our study differs from (llj in two respects: i) 
our Fermi-liquid-like ground state is the result of sponta- 
neous symmetry breaking leading to a "chiral antifcrro- 
magnet" ordering defined in terms of 3-spin products, ii) 
we present analitic results for heat capacity and dynamic 
spin susceptibility in the limit of small Fermi-surface size. 
We consider the Hamiltonian: 
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FIG. 2: (Color online) Flux configurations in Gi,2,3,4 states 
correspondingly. Empty /filled plaquette corresponds to the 
flux +1, — 1 for dodecagons and +i, —i for triangles. 



H= J2 Ji (o-iDi) (o-jn,) + A T P T P ,. (1) 
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Here T p = crf a a^ b cr^ c is the three-spin "exchange" op- 
erator corresponding to the p-th triangle. We take 
Jx,y,z = J, Jx',y',z' = J' and without loss of general- 
ity assume J, J' > 0. Unit vectors n; are parallel to 
x, y and z axis for the corresponding links x, x' , etc. 
Eq. ([1]) reduces to the original Hamiltonian of Ref. [§| 
at A = 0. This spin Hamiltonian is rather special, since 
it posesses large number of independent integrals of mo- 
tion, so called fluxes defined as W p = lX=i CT j, CT j s -i > 
where jo, ji, j n = jo defines a minimal close loop 
on the lattice. All W p commute with Hamiltonian and 
with each other and divide total Hilbert space into sec- 
tors, corresponding to different sets of W p eigenvalues. 
Hamiltonian (fTJ can be diagonalized by means of Ki- 
taev representation of spins via 4 Majorana fcrmions: 
af = icicf, where four Majorana operators Cj, c?, c-', cf 
arc defined on each site of the lattice and satisfy anticom- 
mutation relations | c fi c f| = ^ij^a^- This way, each 
spin (2-dim Hilbert space) is represented by four Majo- 
ranas (4-dim Hilbert space). This is a representation in 
the extended Hilbert space; all physical states should sat- 
isfy the constraint: Di \^ p h ys ) = \^ P hys) for an y lattice 



site i, where Di 



Cicfcfc*. Operators Di is the gauge 



transformation operator for the group Zi . Hamiltonian 
(fTJ), extended to the Hilbert space of Majorana fcrmions, 
reads (iiij = icfcj): 



H 



with t p = 



E '«< 
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Note, that 



uH] = and 



u\j = 1 . These integrals of motion are gauge- dependent; 
they are related to the gauge-invariant fluxes W p = 
n"=i ( — iujsjs-i)- There are two types of W p : 1) fluxes, 
corresponding to the triangular loops W p 3 ^ = ±i, and 2) 



fluxes W, 



(12) 



±1 corresponding to the dodecagon loops. 



They respond differently to the time reversal transfer- 



T< 2 



mation 
shown in Ref 
(fTJ) with A = 

,0) 



w, 



(12) 



but TW, 



(3) 



-Wj, 3) . It was 



12[ that the ground state of Hamiltonian 



corresponds to all W p 



(12) 



T and all 



W p °' are equal (either to i, or to 



-i) ; these two global 



cigenstatcs are related by the T inversion. We show now 
that in some range of couplings J, J' even very small A 
stabilizes another type of the ground state, with vari- 
ables Wp 3 ^ = ±i ordered alternatively (like in the AFM 
Ising model on honeycomb lattice), and with Fermi- line 
of gapless excitations. 

Hamiltonian ([2J can be diagonalized for any periodic 
configuration of the gauge field Uij. However, we restrict 
our consideration to the states with the same flux pe- 
riodicity as the original lattice. Thus we are left with 4 
gaugc-nonequivalcnt states C?i,2,3,4 shown in Fig. [5], plus 
their time- reversal partners GV = TG;. Since the gauge 
field Uij which correspond to the flux configurations G24 
does not fit into 6-site unit cell, in order to describe all 
states G1..4 we use elementary cell containing 12 sites, 
shown in Fig [T] with ei = e x and e2 = e y (hereafter 
lengths are measured in units of the lattice spacing). 

After the gauge is fixed, we are left with the Hamil- 
tonian "H, restricted to the Majorana space, and de- 
noted by H . This Hamiltonian can be diagonalized 
in terms of Fourier-transformed Majorana fields %j) a k = 
, e - c a r , where subscript a = 1..12 enumerates 



I2N 



fcrmionic components inside each of N elementary cells. 
In the Fourier representation the Hamiltonian reads: 



h= E v^k^ 



(3) 



keif 1 



where summation is going over the half of the Bril- 
luen zone K+ = (0 < k x < n, — it < k y < tt). Fourier- 
transformed Majorana fields, restricted to K + , define 
complex fcrmions. Hamiltonian H is diagonal in the 
number of this fcrmions, which results from the trans- 
lational invariance of the system. is a 12 x 12 

gauge-dependent Hermitian matrix. Spectral equation 

dct (ff k - e) = determines twelve bands with disper- 
sions e Qi k- The Fermi-sea energy can then be calculated 
as E = 2 /(ea,k)ea,k, with /(e) being the fermion popu- 
lation numbers (below we imply periodic boundary con- 
ditions). When the ground state of Majorana system in 
some fixed gauge is found, the true ground state of the 
original spin Hamiltonian should be found applying the 
projection operator P — J\ i 1+ 2 P,: . However, if we are in- 
terested in calculation of gauge-invariant quantities (like 
ground state energy or spin susceptibility), there is no 
need for explicit implementation of this projection, and 
calculation can be done in any particular gauge. 

Below we consider vicinity of the point g = g c = s/2 
where the state G2 becomes critical (see below). Ground 
state energies of G1..4 states per unit cell at g = g c are 



at A = 0: E{ 



(0) 



-10.758 J, E { ° ] 



-10.681J, E. 



(0) 



-10.664J, E£ J = -10.610J. That confirms that Gi 
has the lowest global energy at A = 0. However, at 
finite A this energies are simply shifted by ±6A and 
AFM ordcrdcring of t p realized by the G2 state becomes 
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favourable at A > A c = (£ 2 0) - E ( °>) /12 w 6.4 x 10" 3 J. 

The G2 state breaks both T-symmetry and the symme- 
try V of inversion between sublattices of the honeycomb 
lattice, and it can be called "chiral AFM" state. How- 
ever, it is invariant with respect to the combined inver- 
sion "TV. As temperature raises above some critical value 
T c , a phase transition leading to a " chiral-disordered" 
state obeying both T- and P-invcrsions should occur. 
We assume below that T <C T c , and neglect excitations 
which flip chirahties t p . 

Eigenstates of ([3]) are found via matrix diagonaliza- 

tion: Hk = SuH^S^ . Solving the equation detiTk = 
at g = g c we find that zero-energy excitations are located 
at two inequivalent points: Ki and K2. So, there is a sin- 
gle gapless band ei k containing low-energy excitations, 
whereas all other 11 bands £2..i2,k have the gap of the 
order of J. Low-energy excitations are given by 1-k = 
E a S+VvKa+k and <fe,k = Ea <5^a,K 2 +k (hereafter 
for brevity we write S a = S Q .Ki and S a = Si a> jc 2 ). Per- 
turbation expansion up to the second order in k and up 
to the first order in T — g c — g -c 1 leads to the effective 
Hamiltonian of the low-energy excitations: 

H *ff = (^£k02,k - 4>t,k<t>i,k) £k (4) 

|k|«l 

where e k = (3fc^ + k%) - fi , with \i = ^/WJiJT. 

Density of states is defined by J = v J de and is 

equal to v = (2ttJ) 1 . Hamiltonian (jj) determines low- 
energy properties of the spin system ([1]) at small T and 
under condition A > A c . The spectrum is gapful at T < 0. 
Positive r corresponds to the spin metal state. At T <C 
(i < J the heat capacity of spin liquid (per unit cell) 
is C(T) = \T I J , demonstrating standard Fermi-liquid 
behaviour at low temperatures. However, these gapless 
excitations do not carry spin, while spin excitations are 
gapped in this model, as we discuss below. 

We note that similar analysis could be developed for 
the G3 state in the vicinity of the point g = which is 
known to be critical for the G\ state as well However, 
within the model defined by the Hamiltonian (jTJ , the 
state G3 always has energy higher than that of G2 state; 
the latter, however, has large gap near g = \/3, thus we 
work near the point g = \[2 where G2 becomes critical. 

Now we turn to the calculation of frequency-dependent 
spin susceptibility x(cj, T) of the spin metal state. Linear 
susceptibility tensor is proportional to the unit matrix 
due to cubic symmetry of the Hamiltonian ([T]) in the 
spin space. We choose external homogenious field h(t) in 
the z direction which adds the term —h(t) Era °a( r ) to 
the spin Hamiltonian and calculate (a z ). Susceptibility 
reads xH = £r, r ' E Q/ 3 X a pi v ' ~ r > w ) with 

poo 

X*p(T, U )=i e™([** a (r t t),o* fi (p t 0)])dt (5) 
Jo 



according to Kubo formula (r enumerates cells and a, (3 
stay for sites within the same cell, average hereafter 
is taken over the nonperturbed ground state). We 
find, following [?J , that correlation function G a p (r, t) — 

a z a (r , t)a z p (0, 0) ^ is non-zero either for coinciding spins 

or for spins which are connected by z or z' link. This 
means, that G a p(r,t) ~ 6 (r) (since different elemen- 
tary cells are connected by x' , y' links) and in what 
follows wc do not write spatial coordinates explicitly. 
Spin operator creates two Z2 vortices in the neigbour- 
ing plaquettes, which have an excess energy ft a p ~ J, 
thus correlation function G a p(t) oscillates with a fre- 
quency Q a p- Therefore dynamic susceptibility x(w) = 
N Ez (x a/ 3 + Xpp + X a p + Xp a ), where summation goes 
over z and z' links I = (a/3) in the unit cell, contains 
two resonances at frequencies oj — il z , Cl z > . Our goal is 
to find lineshapes of these resonances. 

The sum in parenthesis in the above expression for 
X{u}) can be written as Xi = 4« J °° e iLUt (Gi (t) - G* t (t)) 

with Gi{t) = /e iH Vi e_i ' ff '*V'( + ) 5 where Vz = 

2 (cqj + iu a pcp) is a complex fermion defined on a link 
I. In this expression H[ stays for the Hamiltonian which 
is different from H by inversion of the sign of uf. H[ = 
H + Vi, whereas Vi = 4Ji (tpi~ipi — |). After standard 
algebra we find: 

Gi (t) = ^ (t)T exp (-if Q Vi (r)dr^j ^+ (0)^ . (6) 

The problem of calculation of Gi(t) seems to be simi- 
lar to the Fermi Edge Singularity (FES) problem with 
a separable scatterer. The latter was solved exactly (in 
the infrared limit fit 3> 1) in [l3| by summation of the 
perturbation theory series via the solution of particu- 
lar integral equation. However, our problem is, strictly 
speaking, different, since initial Hamiltonian H is not di- 
agonal in ip^ipi and hence correlation function -Fbz(i) = 
{Ttp l {t)tp l (0)) is not equal to zero identically. However, 
unlike Goi(t) = (T^i (t) xjjf (0)) which has long-time tail 
~ 1/f, the function F ol (t ^> M _1 ) decays very fast with t 
due to exact cancellation between Fermi-surface contri- 
butions coming from different valleys. This means that 
corresponding pairings in the series expansion of (JB]) do 
not lead to any singulary at the threshold and can be ne- 
glected. In this case, the solution is similar to the one pre- 
sented in Ref. [l3| and can be written in terms of the long- 
time (\tfj,\ > 1) asymptotics of Goz —iv + A/<S (i)) 
(here the term with 8(t) is necessary to preserve the cor- 
rect weight J Goi(t)dt, see [l3| for details). The l/t tail 
in the Green function reflects the presense of a band of 
gapless excitations <pi 2 which form Fermi-sea. Reshuf- 
fling of the Fermi-sea by the scattering off the local re- 
pulsive potential Vi leads to the shift of the energy and to 
the change in the power-law exponent in the exact Green 
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FIG. 3: Real (full line) and imaginary (dashed line) parts 
of x(k>',T) at fi/J — 0.05, that corresponds to the exponent 
A = 0.35. 



function, compared to the bare one: 



^ = P e 



(7) 



A-kvcli Ji 



where A/ = 2 (5//7r) - (5j/7r) and 5i = - arctan - 1+4uA Jl ■ 
Frequency 0; ~ J > is the shift in the ground state 
energy due to creation of two fluxes. Obviously, f2; is 
different for z and z' links; finally, £ ~ /j, is the high- 
energy cutoff. Whereas ai is determined by the vicin- 
ity of Fermi-energy only, the parameter A; characterizes 
short-time behaviour of the Green function, and thus is 
determined by the whole spectrum of all 12 fermionic 

bands: a t = £ p , 7 (|S Q7 , P | 2 + |S> 7lP | 2 ) (£ 7lP ) and 



A 



2ii„ 



E P , 7 e 7 , p Im 



iv* z^ Pl7 s 7 ,p— L^p^.pJ ' parameters 
a; and A; are gauge-independent constants, which de- 
pend on the type of the link only. Explicite calcula- 



tion leads to a; 



S a | 



\S, 



\S a \ 



s. 



and 



A; = ci log ( ^ ) +r} u where q = 2u a p Im 



and rii is some number of the order of unity which can 
be found only by numerical integration over K + (it is 
determined by the whole band). Evaluation leads to the 
following result: 



tan <5, = 



-1 



1.04 



0.181n^' 



tan<5 z / = 



1 



0.40 



0.26 In ± 



(8) 

Eq. ((SJ determines phase shifts modulo 7r only. This ambi- 
guity is fixed by the continuity condition: Si(fi — > 0) = 0. 
Since X z < and A 2 ' > for any J/ fx, only \z> diverges 
at the corresponding threshold (while x z still has a cusp) ; 
therefore below we concentrate on the contribution of the 
z' links only, \ z i = A. Using the results (|7|8p . we calcu- 
late spin susceptibility close to the fi z / resonanse, with 
uj' = u - SI,,. We find Xt=o (<«/) = A o r (A) (-Co/w') A , 



where Xo 



AN 
3tt 



J . Note, that = below the thresh- 



old (u/ < 0) as it should be at T = 0. These results can 



be easily generalized to the finite temperature T -C fi- 
As was shown in [l4j , finite-temperature correlation func- 
tion in the FES problem can be obtained from the zero- 
temperature one by substitution t — > sin ^ rt . For the 
susceptibility, that gives: 



X ) = Xo 



\2irT J r(i±A-^) 



(9) 



This function is plotted in Fig. [3] The ma- 
jor effect of nonzero temperature is the appearence 
of absorption below threshold: xi^^'^T) « 

Xt=o ( w ') fl + isinvrAe - '" '/ T ^ ; in addition, the reso- 
nant peak appears to be smeared out: xd^'l < T) s; 

v ( Co A 

Xo e V 27rT / ' 

In conclusions, we have shown that certain (numer- 
ically, very weak) modification of the Yao-Kivelson ver- 
sion of Kitaev spin lattice leads to the ground-state of the 
Fermi-liquid type, with a Fermi energy /x cx v2 — J' /J. 
We have studied the model in the continuum limit of 
small Fermi-circle \i -C J and at low temperatures T <C 
fx. Gapless excitations of the Fermi-sea do not carry spin 
themselves, but they determine the shape ((9j of the res- 
onance peak in the dynamic spin susceptibility. 
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